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Abstract [Background]: The quantification of uncertainty has become a common requirement in reactor physics 


analysis, with the covariance data of nuclear data serving as the foundational data for conducting uncertainty 


quantification. 


[Purpose]: This study aims to develop a covariance data generation module, covar calc, in the nuclear 


data processing software NECP-Atlas to produce continuous energy covariance data for Monte Carlo programs and 


multi-group covariance data for deterministic programs. [Methods]: Covar calc module was developed to process all 


covariance data provided in evaluated nuclear databases, according to the different storage formats of nuclear data and 


different computational methods. This included the covariance data of various parameters such as average fission 
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neutron multiplicities, cross sections, angular distributions of secondary particles, fission spectrum, resonance 
parameters, and neutron activation cross-sections. [Results]: A comparative verification is carried out with the 
covariance data production module in the nuclear data processing software NJOY21 to demonstrate the accuracy of the 
multi-group covariance databases produced by NECP-Atlas and the maximum bias is less than 0.196. Utilizing Monte 
Carlo calculation code NECP-MCX and uncertainty analysis code NECP-UNICORN, sensitivity coefficients for 
different benchmarks are calculated. By incorporating both continuous energy covariance databases and multi-group 
covariance data, and utilizing the "Sandwich formula," the final uncertainties are computed. The results indicate that the 
uncertainties calculated using the multi-group covariance data generated by NECP-Atlas exhibit comparable accuracy 
to those obtained with NJOY21. [Conclusions]: The precision in creating both continuous energy covariance and multi- 
group covariance meets the requirements for usage in Monte Carlo programs and deterministic programs, validating the 


efficacy of covar calc module within NECP-Atlas for uncertainty quantification in reactor physics analysis. 


Key words Covariance data, Uncertainty quantization, Evaluated nuclear data, Nuclear data process, Covariance 


calculation 


随 着 高 保 真 核反应 堆 物 理 计 算 方 法 的 发 展 ， 由 于 理 六 从 模型 等 因素 带 来 的 不 确定 性 越 来 越 小 。 核 数据 作为 核 
反应 堆 物理 计算 最 基础 的 输入 参数 ， 其 理论 计算 、 实 验 测 量 、 数 据 评价 等 过 程 ， 不 可 避免 的 产生 一 定 的 不 
确定 度 吊 ， 其 不 确定 度 已 经 成 为 核反应 堆 物 理 计算 结果 不 确定 度 的 重要 来 源 , 因此, 进行 核反应 堆 高 精度 模 
拟 给 出 最 佳 估 计 值 并 量化 核 数据 造成 的 计算 结果 的 不 确定 度 已 成 为 核反应 堆 设 计 的 普遍 需求 。 

核 数据 的 不 确定 性 以 协 方差 的 形式 保存 在 评价 和 数据 库 中 。 目 前 ， 评 价 核 数据 库 给 出 了 平均 裂变 中 子 
数 、 共 振 参 数 、 中 子 反 应 截面 、 次 级 粒子 角度 分 布 、 裂 变 谱 、 活 化 截面 等 数据 的 协 方差 中。 随 着 核 数 据 不 确 
定 度量 化 需求 的 不 断 增加 ， 评 价 核 数 据 库 中 提供 的 协 方差 数据 也 在 大 幅 增 加 。 美 国 的 ENDF/B 系列 评价 核 
数据 库 ， 从 ENDF/B-VILOP!, ENDF/B-VIL1U5] ENDF-B/VIIL0G5]， 提 供 协 方差 数据 的 核 素 数量 以 及 协 方差 
数据 的 类 型 均 有 了 明显 的 增加 ， 目 前 最 新 的 ENDF/B-VIIL0 评价 核 数据 库 提供 了 182 个 核 素 的 评价 协 方差 
数据 库 ， 中 国 核 数据 中 心 最 新 发 布 的 CENDL-3.2 中 ， 提 供 了 77 个 核 素 的 中 子 截面 协 方差 数据 器 
评价 核 数据 库 提 供 的 协 方差 数据 〈 下 文 称 为 评价 协 方差 数据 ) 具有 类 型 多 样 、 反 应 道 众 多 、 协 方差 网 
格 不 统一 、 相 互 关系 复杂 等 特点 ， 不 能 直接 应 用 于 核反应 堆 物 理 计算 的 不 确定 度量 化 ， 因 此 需要 先 对 评价 
核 数 据 库 中 的 协 方差 数据 进行 处 理 ， 制 作为 特定 格式 协 方差 数据 库 。 目 前 国际 上 能 够 基于 对 评价 协 方差 数 
据 进 行 处 理 的 程序 较 少 。 国 际 上 广泛 使 用 的 核 数据 处 理 程序 NJOY 在 2007 年 接 入 了 日 本 原子 能 机 构 开 发 
的 ERRORJ 程序 中 ， 形 成 了 协 方差 数据 产生 模块 ERRORR， 可 计算 多 群 协 方差 ， 但 其 对 次 级 粒子 角度 分 布 
的 处 理 仅 限于 弹性 散射 反应 道 的 一 阶 勒 让 德 系数 外。 国内 的 核 数据 处 理 程 序 如 中 国 核 数 中 心 开 发 的 RULER 
程序 名、 清华 大 学 开发 的 RXSP 程序 00 等 不 具备 处 理 评 价 协 方差 数据 的 功能 。 从 实现 我 国 核电 软件 自主 化 
的 角度 而 言 ， 研 发 自主 化 的 协 方差 数据 计算 模块 的 功能 是 十 分 必要 的 。 男 一 方面 ， 使 用 蒙特 卡 罗 程 序 计 算 
连续 能 量 灵敏 度 系数 的 理论 与 功能 已 经 得 到 了 逐步 开发 并 完善 04.2， 蒙 特 卡 罗 程 序 因 为 几何 适应 性 强 ， 计 
算 精度 高 的 优点 ， 在 复杂 几何 的 条 件 下 计算 得 到 的 复杂 响应 的 灵敏 度 系数 与 传统 确定 论 方法 的 计算 结果 有 
不 同 程 度 上 的 区 别 L3]。 目 前 ， 国 际 上 核 数 据 处 理 程序 缺乏 产生 适用 于 蒙特 卡 罗 程 序 的 连续 能 量 协 方差 数据 
的 产生 功能 ， 因 此 有 必要 开发 产生 连续 能 量 协 方差 的 功能 ， 满 足 基 于 蒙特 卡 罗 计 算 程 序 的 不 确定 度量 化 需 
求 。 


本 文 在 西安 交通 大 学 核 工程 计算 物理 (Nuclear Engineering Computational Physics, NECP) 实验 室 自主 研 
发 的 核 数据 处 理 程序 NECP-Atlas 的 基础 上 04, 开发 了 应 用 协 方差 数据 产生 模块 , 实现 了 利用 评价 协 方差 数 
据 计 算 生成 连续 能 量 格 式 协 方差 数据 和 多 群 格式 协 方差 数据 (以 下 将 两 者 统称 为 应 用 协 方 差 数 据 ) 的 功能 ， 
并 对 该 模块 进行 了 验证 。 
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1 ”应 用 协 方差 数据 计算 方法 

目前 , ENDF-6 格式 评价 核 数 据 库 中 给 ed ela ti MM 
出 射 粒子 能 量 分 布 、 中 子 活化 截面 的 协 方差 数据 外 。 以 下 分 别 介绍 其 对 应 的 应 用 协 方差 数据 计算 方法 。 

11 平均 裂变 中 子 数 、 中 子 反 应 截面 、 中 子 活化 截面 应 用 协 方 差 的 计算 


平均 裂变 中 子 数 、 中 子 反 应 
本 文 以 中 子 反 应 截 


SA 


fi PL 4 — EE, Z8 T Je RA A 235 33008 A AAE, 
反应 道 的 能 量 网 格 进行 合并 ， 构 建 统一 的 联合 能 量 网 格 。 此 外 ， 评 价 协 方差 数据 的 
要 输出 的 多 群 能 群 结构 边界 点 通常 也 是 不 一 致 的 ， 因 此 还 需要 将 用 户 输入 的 多 群 角 
合 能 量 网 格 。 联 合 能 量 网 格 的 构建 过 程 如 图 1 所 示 。 
Emin Ei Emax Emin E2 E3 Emax Emin E, 
E, E, 
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E; E; 
E 3 E 3 
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Energy grid of reaction 1 Energy grid of reaction 2 
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H. 


in MF 33 in MF 33 


过 数学 关系 


一 致 ， 


E 导 得 到 。 


平价 协 方差 数据 有 两 种 格式 , 分 别 为 NI 格式 (NI-type) 和 NC 格式 (NC-type)。NI 格 式 为 显 式 格式 , 评价 
核 数 据 库 中 直接 给 出 当前 反应 道 与 其 他 反应 道 之 间 的 协 方差 数据 ; NC 格式 为 隐 式 格式 , 评价 核 数 据 库 未 直 
接 给 出 协 方差 数据 ， 当 前 反应 道 的 协 方差 根据 分 反应 道 的 协 方差 数据 通 


Sk E. 


用 5R 


评价 协 方差 数据 通常 分 为 若干 能 量 段 进行 存储 ， 各 能 量 段 协 方差 数据 不 同 ， 并 且 不 同 反 应 道 的 协 方差 
需要 将 评价 核 数 据 库 中 的 所 有 
网 格 和 用 户 最 后 需 
E 群 结构 的 边界 点 并 入 联 


E; E3 Emax 


Union energy grid 


图 1 基于 两 个 不 同 反应 道 协 方差 能 量 网 格 的 联合 能 量 网 格 构建 过 程 


Fig. 1 Construction process of union energy grid based on the grids of two different reactions- 
联合 能 量 网 格 上 的 通 量 和 截面 通过 下 式 计算 : 
É, v 人 (已 Ja (1) 
NOLO é(E)x(E)dE zb ((E)x(E)dE 
S MEINES (2) 


uc T 


7 


式 中 : 7 为 联合 能 量 网 格 序号 ，# 为 第 7 网 格 的 通 量 ，x, 为 第 7 网 格 的 截面 ，G(E) 为 能 量 点 处 的 通 量 ; 
x(E) 为 能 量 点 5 处 的 截面 值 ， 由 NECP-Atlas 经 共振 重 构 、 多 普 勒 展 宽 等 计算 等 模块 产生 0519。 其 中 yp(E) 
可 以 由 用 户 选择 , 也 可 以 使 用 NECP-Atlas 计算 得 到 的 通 量 .将 式 (2) 离 散 化 ， nmi iot P oa, 
I 
得 到 : 
um 
el -~ 27 OX : (3) 
iel 
HÀ mS mm \ 式 计算 。 
由 于 NC 格式 协 方差 数据 是 由 NI 格式 协 方差 数据 线性 组 合 得 到 ， 因 此 需要 首先 在 联合 网 格 上 计算 NI 
格式 的 协 方差 数据 。 对 于 某 反 应 道 在 某 能 量 点 处 的 截面 x 和 男 一 反应 道 在 男 一 能 量 点 处 的 截面 y, 评价 核 数 
据 库 中 可 能 给 出 多 个 独立 的 NI 格式 协 方差 数据 ， 因 
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nim 


cov( X, y) J= Zost (x, y), (4) 


式 中 : nim 为 当前 两 反应 道 评 价 协 方差 中 包含 的 NI 格式 的 个 数 。 cov (Xx,y) 为 第 n 个 NI 格式 协 方差 数据 ， 
根据 式 (3) 和 式 (4) 以 及 协 方差 数据 的 传递 性 质 071， 在 联合 能 量 网 格 上 ， 某 有 反应 道 第 I 个 网 格 下 的 截面 与 另 
一 反应 道 第 J 个 网 格 下 的 截面 yj 的 协 方差 可 以 写 为 : 


nim nim 


cov(x,, y,) Daa CERED y cov (xy. Yy), (5) 


n=l iel 
jeJ 


RP: Xp Yy 分别 为 某 反 应 道 第 1 网 格 内 第 i 个 能 量 点 处 的 截面 值 ， 以 及 另 一 反应 道 第 J 网 格 内 第 7 个 能 量 
点 处 的 截面 值 。cov (xj, yy) 为 截面 x 和 截面 ， 在 第 个 NI 格式 中 的 协 方差。 
在 联合 能 量 网 格 上 ，%, 和 的 相对 协 方差 计算 公式 为 : 


nim 


uut P 5 
Xy, Xy; 

在 上 述 联合 网 格 下 的 协 方差 数据 的 计算 过 程 中 ， 若 在 式 (1) 中 选择 中 子 通 量 选择 常数 谱 ， 此 时 计算 得 到 
的 协 方差 数据 与 问题 无 关 ， 获 得 的 协 方差 数据 库 即 为 连续 能 量 协 方差 数据 。 由 于 联合 能 量 网 格 内 所 有 能 量 
点 的 协 方差 相同 ， 因 此 连续 能 量 协 方差 数据 库 中 对 某 一 能 量 网 格 仅 给 出 一 个 协 方差 数据 。 当 蒙特 卡 罗 程 序 
需要 使 em UE caus 直接 使 用 该 能 量 点 所 在 网 格 的 协 方差 数据 。 

若 需 要 得 到 特定 能 群 结构 的 多 群 协 方差 数据 ， 还 需要 将 联合 能 量 网 格 上 的 协 方差 数据 按 多 群 能 群 结构 
进行 归并 。 群 通 量 和 多 群 截面 计算 公式 为 : 


gl 


rcov (x, y, 


fk ES > 力 (7) 
IeK 
2 x, 
Xy — IeK (8) 


K 
式 中 : 从 为 多 群 结构 中 上 第 KK 群 的 通 量 ，Y) 为 联合 网 格 上 第 了 个 网 格 的 通 量 ，Xk 为 多 群 结构 中 某 反 应 道 
区 第 天 群 的 多 群 截面 值 。 同 样 ， 另 一 反应 道 了 的 第 过 群 多 群 截面 可 以 写 为 ; 


27 


了 = yeL (9) 
$, 
根据 协 方差 的 可 加 性 ，Xx MY, 之 间 的 多 群 协 方差 可 以 由 联合 网 格 上 的 协 方差 归并 得 到 : 
> pp, cov (x,y) 


cov(X p Y,)= EET C (10) 
KTL 


将 以 上 绝对 协 方差 除 以 多 群 截面 值 后 得 到 相对 协 方差 : 
P3 pp, cov(x,. y,) 


rcov ( X ,.. Y, EM 11 
( 中 $6, X kY, e 


对 于 NC 格式 协 方差 ， 评 价 核 数据 库 中 给 出 了 目标 反应 道 与 其 他 反应 道 之 间 的 线性 关系 ， 例 如 ， 菜 有 反 
应 道 第 了 网 格 和 另 一 反应 道 第 /网 格 上 的 截面 值 和 ， 可 以 由 其 他 反应 道 通过 如 下 线性 加 和 获得 : 
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= $C(a)A, (a) (12) 
y, - C(b)B, (b) (13) 


RP: C(a) 为 目标 反应 道 和 某 一 反应 道 a 之 间 的 线性 系数 ，Ai (a) 为 反应 道 a 在 联合 能 量 网 格 上 第 7 网 
KRH: C(b) 5 B, (b) 3x 3.5 C(a) RV A, (a) 类 似 。 将 式 (12) 和 (13) 代 入 式 (11), 利用 协 方差 数据 的 
传递 关系 ， 得 到 多 群 结构 上 NC 格式 相对 协 方差 数据 : 
2 > C,( b)6,ó, cov| A, (a). B, (P) | 
cov( X ,了 i (14) 
( T 中 9.0, X kY, 
终 ， 在 输出 的 多 群 能 群 结构 中 ， 多 群 协 方差 为 NI 格式 中 与 NC 格式 中 计算 得 到 的 协 方差 之 和 。 多 群 
XE pA. 多 群 活化 截面 的 协 方差 数据 同 理 计算 得 到 。 


1.2 次 级 粒子 角度 分 布 协 方差 数据 计算 


通过 处 理 次 级 粒子 角度 分 布 的 协 方差 数据 ， 计 算得 到 多 群 勒 让 德 系数 的 协 方差 数据 和 平均 散 财 角 余弦 
的 协 方差 数据 。 多 群 勒 让 德 系数 的 协 方差 数据 与 中 子 截面 的 处 理 流程 相似 ， 同 样 需要 先 在 联合 网 格 下 计算 
协 方差 ， 再 以 核反应 率 为 权重 ,将 协 方差 归并 至 输出 的 多 群 能 群 结构 上 ， 本 节 不 再 次 述 。 以 下 推导 多 群 平 
均 散射 角 余 弦 协 方差 的 计算 公式 : 

实验 室 系 下 平均 散射 角 余 弦 的 计算 公式 可 以 写 为 : 


m =f uf (pn) dn (15) 


RP: M 为 实验 室 系 下 平均 散射 角 余 弦 ， Lu, 为 实验 室 系 下 散射 角 余弦 f (uu) THO SER IU ORE 
的 概率 分 布 。 
将 了 (4) 用 勒 让 德 多 项 式 展开 : 


f (1)=>, 2 a, (E)P,(u) (16) 


RP: NL 为 勒 让 德 多 项 式 系数 个 数 ，w , (E) 为 实验 宝 系 下 勒 让 德 多 项 式 系数 。 将 式 (16) 代 入 式 (15)， 同 
时 利用 勒 让 德 多 项 式 在 [-1,1] 上 的 正 交 性 09， 式 (15) 可 以 表示 为 : 


— p X 2n+1 
=) P(Y 28a, (E)P, (n) dmi =a, (E) an 


n=0 


上 式 表 明 ， 实 验 室 系 的 平均 散射 角 余 弱 与 实验 室 系 的 一 阶 勒 让 德 系数 相等 ， 因 此 其 协 方差 也 相等 ， 同 样 地 
在 联合 网 格 下 协 方差 也 相等 ， 即 : 


cov (zB) = cov(a,,,.4,,., ) (18) 


Xu. A. HQQSSURSUR FREA HERR TIRARA: a a 为 实验 室 系 下 不 同 入 射 能 
的 一 阶 勒 让 德 系数 。 当 评价 核 数据 库 中 直接 给 出 实验 室 系 下 的 勒 让 德 系数 时 , 其 协 方差 数据 可 以 直接 使 用 ， 
当 评价 核 数据 库 中 给 出 质心 系 下 的 勒 让 德 系数 时 ， 其 协 方差 数据 还 需要 进行 相应 的 坐标 变换 ， 得 到 实验 室 
系 下 的 勒 让 德 系数 09, 最 终 可 以 推导 出 实验 室 系 下 平均 散射 角 余 弦 和 质心 系 下 的 勒 让 德 系数 之 间 的 关系 为 ; 
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| zt | Qr , a9 
MZR SR2 ) 和 27 -1 2n+3\R R i: 
` ` Ah L^. ^L B NA JN e (A«1)(-Q) ` 
WO ME MN UA 


量 , O 为 质心 系 下 的 反应 阔 能 , 已 为 质心 系 下 的 入 射 中 子 能 能 量 . 将 式 (19) 中 心 ,前 的 系数 , 例如 、 I- 
等 定义 为 ， 则 在 联合 能 量 网 格 下 , 由 质心 系 下 勒 让 德 系数 计算 实验 室 系 下 平均 散射 角 余弦 的 协 方差 的 公 


式 为 : 


X 


cov( 7 ,, d, )= Se cov(a d s ajy) (20) 


i,j=1 
AP: aur 0,5, 分 别 为 质心 系 下 的 第 1 网 格 内 第 i 阶 勒 让 德 系数 与 第 J 网 格 内 第 j 阶 勒 让 德 系数 。 
式 (19) 和 (20) 计 算 过 程 中 需要 用 到 联合 网 格 下 的 勒 让 德 系数 , 联合 网 格 下 的 勒 让 德 系数 计算 公式 与 式 (2) 


EIE a(E)dE BU (E)dE 
je(E)dE fh 


(21) 


AP: aj 为 第 了 网 格 内 的 某 一 阶 数 下 的 勒 让 德 系数 。 式 (21) 需 要 用 由 式 (1) 计 算得 到 的 中 子 通 量 。 在 中 子 通 
量 计算 时 选择 常数 详 ， 此 时 计算 得 到 的 平均 散射 角 余弦 协 方差 数据 与 问题 无 关 ， 输 出 的 协 方差 数据 库 即 为 
连续 能 量 协 方差 数据 ， 可 直接 用 于 蒙特 卡 罗 程 序 的 不 确定 度量 化 。 

对 于 多 群 平均 散射 角 余弦 协 方差 的 计算 ， 以 核反应 率 为 权重 ,将 归并 至 输出 网 格 上 ， 类 似 于 式 (11), F 
均 散射 角 余弦 相对 协 方差 计算 公式 如 下 : 


3 


2" $,0,0,0, cov( Z; Hg ) 


IeK 
JeL 


rcov( 7 y. hy)= Q2) 


$,0,0,0, «HL 


式 中 :五 ， 厂 ,一 一 多 群 能 群 结构 上 第 大群 和 第 工 群 平均 散射 角 余 蓄 ，ar 0, 一 一 多 群 能 群 结构 上 第 天 
群 和 第 工 群 对 应 反应 道 的 中 子 截面 0, . 0, 一 一 联合 网 格 上 的 第 了 个 网 格 和 第 个 网 格 对 应 反应 道 的 中 子 


13 裂变 谱 多 群 协 方差 数据 计算 
评价 核 数据 库 给 出 了 裂变 反应 道 出 射 粒 子 能 量 分布 的 协 方差 ， 即 裂变 谱 的 协 方差 数据 。 对 于 裂变 谱 协 
方差 数据 的 计算 ， 首 先 在 联合 网 格 上 计算 裂变 谱 的 协 方差 数据 ， 然 后 以 中 子 通 量 为 权重 ， 将 协 方差 数据 归 
并 至 输出 的 能 群 结构 上 ， 与 式 (10) 类 似 。 当 入 射 中 子 能 量 为 E 时， 裂变 谱 协 方差 数据 的 计算 公式 如 下 : 

2 Ahio E Ei) p (E-E) 


cov( p(E > Ez). p(E E; )- fe T (23) 


: p(E E,), p(E 一 ) 一 一 入 射 能 量 为 E 的 中 子 ， 散射 至 输出 能 群 结构 第 KK 群 和 第 工 群 的 概 


X 


Z; p(E— E), p(E — EL) ——P fits E 的 中 子 ， 散射 至 联合 网 格 上 第 了 个 网 格 和 第 / 个 网 格 的 
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概率 。 


在 计算 多 群 裂 变 谱 的 相对 协 方差 数据 之 前 ， 需 要 计算 入 射 能 量 为 妃 时 的 多 群 裂 变 谱 ， 计 算 公 式 如 下 : 
2,0 esu 
yA A 
2,4, 2,0, 
8 8 


式 中 ，8' 、g 一 分 别 为 入 射 能 群 和 出 射 能 群 ，。 — 第 g 群 多 群 裂变 谱 ，a/ ，， — ERRIRE 


ME: el. eh 入 射 能 量 的 下 限 和 上 限 。 
与 截面 的 相对 协 方差 计算 类 似 ， 在 裂变 谱 多 群 协 方差 数据 基础 上 除 以 多 群 裂变 谱 ， 得 到 裂变 谱 相 对 协 
方差 的 计算 公式 如 下 : 


,8' el(el,eh) (24) 


>. 44,cov(p(E > Ej). p(E > E;)) 


IcK 
rcov( p(E > E; ), p(E > E;))- = (25) 
KTLAK^AL 


14 ”共振 参数 多 群 协 方差 数据 计算 


评价 核 数 据 库 给 出 了 共振 参数 之 间 的 协 方差 数据 ， 需 要 将 其 转化 为 共振 反应 截面 之 间 的 协 方差 数据 ， 
并 与 中 子 截面 协 方差 数据 加 和 获得 共振 能 区 截面 的 协 方差 ， 用 于 不 确定 度量 化 。 与 处 理 中 子 截面 协 方差 数 
据 过 程 不 同 ， 共振 参数 协 方差 数据 处 理 时 ， 直接 将 共振 参数 协 方差 转化 为 输出 网 格 上 共振 反应 道 的 多 群 截 
面 的 协 方差 。 转 化 公式 如 下 : 
JE Jew, ) (26) 
D, 


cov(X,, X,)— XE: 


RP: 人 J—— REKREA X, X,— 一 不 同 能 群 的 共振 反应 道 的 截面 ，T, 、T ,一 一 


不 同 的 共振 参数 ， 评价 核 数 据 库 中 提供 了 不 同 共振 参数 之 间 的 协 方差 cov, D) ， 对 于 灵敏 度 系数 
OXj /ol';， 采 用 直接 扰动 方法 进行 计算 。 

在 可 分 辨 共振 能 区 ， 评 价 核 数 据 库 对 于 不 同 的 共振 公式 ， 分 别 给 出 了 共振 峰 能 量 、 不 同 反 应 道 共振 峰 
宽度 和 散射 半径 的 协 方差 数据 。 其 中 截面 对 共振 峰 能 量 和 不 同 反应 道 共 振 峰 宽度 的 灵敏 度 系数 ， 直 接 扰动 
法 计算 公式 如 下 : 

(e E "(E)-co (E 
OT,  Ag-te GT, Age 2AT, 
RP: ar (E). o^ (E) — 在 能 量 巨 处 对 第 ;个 共振 参数 扰动 T, AT, RT; - AT, Je tf ARTo. 


不 同 共振 公式 扰动 的 共振 参数 不 同 ， 如 对 于 Single-Level Breit- -Wigner 公式 和 Multilevel Breit-Wigner 公式 ， 
分 别 扰动 共振 峰 能 量 、 反应 道 宽度 、 俘 获 反应 道 宽 度 和 裂变 反应 道 宽 度 、 散 射 半径 ， 对 于 Reich- 
Moore 公式 ， 分 别 扰动 共振 峰 能 量 、 弹 性 散射 反应 道 宽度 、 俘 获 反应 道 宽 度 、 第 一 部 分 裂变 反应 道 宽 度 和 
第 二 部 分 裂变 反应 道 和 径 。 其 中 ， 对 共振 峰 能 量 扰动 0.01%， 对 其 他 核 数 据 扰动 196. 


2 ”程序 开发 及 应 用 协 方差 数据 库 直 接 数 值 验证 


基于 上 述 理论 ， 在 NECP-Atlas 开发 了 协 方 差 数 据 产生 模块 covar cale, 实现 了 基于 评价 协 方差 数据 加 
工 制作 应 用 协 方差 数据 的 功能 。 以 截面 的 应 用 协 方差 数据 为 例 ， 其 产生 过 程 如 图 2 所 示 。NECP-Atlas 首先 
读 取 评价 核 数 据 及 用 户 输 入 ; 然后 通过 共振 重 构 、 多 普 勒 展 宽 等 模块 ， 加 工 生 成 特定 温度 点 下 的 点 截面 数 


sk 
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据 库 (PENDF); 根据 评价 协 方差 数据 ， 构 建 联 合 能 量 网 格 ， 如 果 最 终 需 要 输出 多 群 协 方差 数据 库 ， 同 时 将 
多 群 能 群 结构 边界 点 并 入 联合 能 量 网 格 ， 基 于 联合 能 量 网 格 ， 计 算 各 网 格 上 的 截面 、 通 量 、 协 方差 数据 ; 
最 终 ， 对 于 连续 能 量 协 方差 数据 ， 直 接 输 出 计算 结果 ， 对 于 多 群 协 方差 数据 ， 则 需要 将 联合 能 量 网 格 下 的 
结果 归并 至 多 群 能 群 结构 上 。 


Read Evaluated nuclear data 
file 


Read User's input 


Reconstruct, doppler 
broadening module, generate 
PENDF 


Construct Union energy grid 
and construct NC Coefficient 
matrix 


Merge multi-group 
boundary points into 
union energy grid 


enerate Conti nus 
energy covariance 2 


Yes 


Calculate cross sections, flux, 
covariance data in union energy 


grid 


: Merge cross sections 
enerate Continus No d Bi 
: flux, covariance data into 
energy covariance 2 : 
multi-group structure 
Yes 


Output multi-group 
covariance data 


Output continuous energy 
covariance data 


图 2 应 用 协 方差 数据 产生 过 程 


Fig.2 The process of applying covariance data generation 


基于 covar cale, 利用 ENDF/B-VIILO 评价 核 数据 库 ， 产生 了 相应 的 多 群 协 方差 数据 库 和 连续 能 量 协 方 
差 数 据 库 。 为 验证 covar_calc 产生 的 应 用 协 方差 数据 库 的 正确 性 ， 本 文 将 NECP-Atlas 产生 的 多 群 协 方差 数 
据 库 与 核 数 据 处 理 程序 NIOY21 在 相同 条 件 下 产生 的 协 方差 数据 库 进 行 直接 对 比 ， 构 造 相对 偏差 矩阵 ， 统 
计 各 相对 偏差 矩阵 中 元 素 的 最 大 值 和 矩阵 的 F 范 数 〈 和 矩阵 中 所 有 元 素 的 平方 和 的 开 方 值 )， 通过 这 两 个 数值 
的 大 小 来 反映 多 群 协 方差 数据 库 的 偏差 程度 。 
在 计算 过 程 中 ， 共 振 参 数 的 协 方差 数据 通过 转化 得 到 中 子 截面 的 协 方差 数据 。 因 此 通过 直接 对 比 中 子 
截面 的 多 群 协 方差 数据 库 ， 可 以 同时 验证 NECP-Atlas 处 理 共振 参数 协 方差 和 中 子 截 面 协 方差 数据 的 功能 。 
表 1 给 出 了 本 文 计 算得 到 的 不 同 核 素 、 不 同类 型 核 数 据 的 相对 偏差 矩阵 的 最 大 值 和 下 范 数 ， 对 平均 裂变 中 
子 数 、 中 子 截 面 〈 含 共振 参数 )、 平 均 散 射 角 余 弦 、 裂 变 谱 等 不 同类 型 核 数 据 进行 了 验证 。 
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表 1 不 同类 型 核 数 据 的 协 方差 数值 验证 结果 


Tablel Numerical verification results of the covariance of different nuclides and different nuclear 


data types 
核 素 核 数 据 1 核 数 据 2 最 大 值 F 范 数 核 素 核 数据 1 BOR 最 大 值 F 范 数 
gii Nuclear Nuclear Maximum Frobenius Wi Nuclear Nuclear Maximum Frobenius 
Nuclides Nuclides 
data 1 data 1 value /% norm data 1 data 1 value /% norm 
-]1.423E- 7.315E- -1.226E- 5.150E- 
U-235 Oclas Oclas U-235 Oy Oy 
04 07 01 05 
加 - -6.145E- 5.054E- - E -6.674E- 5.102E- 
U-235 Vi V U-235 Vi V 
05 07 $ 05 07 
n - 2.373E- 2.398E- -2.683E- 6.779E- 
Pu-239 n H U-235 X X 
02 05 04 05 
= - -3.728E- 5.414E- -1.973E- 7.885E- 
O-16 a u U-238 X X 
02 05 04 05 


从 上 表 可 以 看 出 ， 对 于 所 测试 的 反应 道 ， NECP-Atlas 与 NJOY21 的 相对 偏差 矩阵 中 的 最 大 偏差 范围 为 
-0.1226%~0.02373%， 相 对 偏差 矩阵 的 下 范 数 不 超过 8X105， 证 明 NECP-Atlas 产生 的 协 方差 数据 的 正确 


性 。 


3 ， 协 方差 数据 库 在 不 确定 度量 化 中 的 验证 


为 进一步 验证 多 群 协 方差 数据 计算 的 正确 性 ， 分 别 基于 NECP-Atlas 和 NJOY21 产生 了 69 群 多 群 协 方 
差 数据 库 ， 采 用 NECP-UNICORNP9o 程 序 基 于 以 上 多 群 数据 库 库 计算 了 UAM-LWRs 系列 基准 题 中 的 TMI- 
1 燃料 栅 元 基准 题 P20， 量 化 了 不 同 核 素 各 反应 道 对 无 限 增殖 因子 kir 造成 的 相对 不 确定 度 ， 其 计算 公式 见 
(28): 


Alii 100% (28) 


inf 
导 到 的 无 限 增殖 因子 kinr 的 不 确定 度 。 


d 2 给 出 了 热 态 满 功率 工 况 下 的 TMI-1 栅 元 中 不 同 核 素 不 同 反 应 道 对 kine 造成 的 的 相对 不 确定 度 ， 对 
比 了 基于 NECP-Atlas 与 NJOY21 产生 的 多 群 协 方差 数据 库 的 计算 结果 的 偏差 , 从 表 中 可 以 看 出 NECP-Atlas 
产生 的 多 群 协 方差 数据 库 ， 具 有 与 NJOY21 相当 的 精度 ， 不 确定 度 的 相对 偏差 全 部 在 0.1% 以 下 。 
表 2 基于 多 群 协 方差 的 TMI-1 栅 元 无 限 增殖 因子 kw 相对 不 确定 度 计算 结果 


Table2 Calculation results of TMI-1 cell relative uncertainties of infinite multiplication factor 
based on multi-group covariance 


式 中 : A at 一 一 NECP-UNICORN it 4 


ENS 


核 素 及 反应 ”NJOY21 结 NECP-Atlas A 核 素 及 反应 ”NJOY21 结 — NECP-Atlas i 
. 相对 偏差 /% 相对 偏差 /% 
道 果 /% 结果 /9% 道 果 /% 结果 /% 
realative realative 
Nuclides and NJOY21 NECP-Atlas Nuclides and NJOY21 NECP-Atlas 
bias/96 bias/96 
recations results/26 results/% recations results/% results/% 
U-235 Gelas 1.172E-03 1.172E-03 0.000 U-238 Gelas 1.410E-01 1.410E-01 -0.001 
U-235 ainel 1.454E-03 1.454E-03 0.000 U-238 ainel 1.253E-01 1.253E-01 0.000 
U-235 on 2.715E-04 2.715E-04 0.004 U-238 o». 1.198E-02 1.199E-02 0.084 
U-235 ot 7.638E-02 7.638E-02 0.000 U-238 ot 1.564E-02 1.564E-02 0.000 
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U-235 o, 1.967E-01 1.967E-01 0.000 U-238 o, 2.906E-01 2.906E-01 0.001 
U-235 v 1.251E-01 1.251E-01 0.000 U-238v 71.336E-02 7.336E-02 0.000 
H-1 Gelas 3.278E-02 3.278E-02 0.000 O-16 Gelas 1.996E-02 1.996E-02 0.000 
H-1 oy 9.495E-02 9.495E-02 0.000 O-16 oy 1.998E-04 1.998E-04 0.000 
为 验证 NECP-Atlas 中 连续 能 量 协 方差 数据 制作 功能 的 正确 性 ,使 用 蒙特 卡 罗 程 序 NECP-MCXP223]， 计 
算 了 ICSBEP 基准 题库 中 的 Godiva 基准 题 、Jezebel 基准 题 P4 中 主要 核 素 对 有 效 增 殖 因 子 ker 的 灵敏 度 系 
数 ， 并 基于 连续 能 量 协 方差 数据 计算 得 到 相对 不 确定 度 结果 ， 与 NECP-UNICORN 基于 多 群 协 方差 数据 计 
算 的 结果 进行 了 对 比 。 在 基准 题 Godiva 基准 题 及 Jezebel 基准 题 中 ， 对 比 核 素 分 别 选 为 U-235 及 Pu-239, 
不 确定 度 计算 结果 如 表 3 所 示 。 从 表 中 可 以 看 出 ， 基 于 连续 能 量 协 方差 数据 的 相对 不 确定 度 计算 结果 和 多 
群 协 方差 计算 结果 最 大 偏差 为 1.45%， 但 此 时 绝对 偏差 仅 为 4X105 (4 pcm)， 证 明了 NECP-Atlas 中 连续 
能 量 协 方差 制作 功能 的 正确 性 ， 可 以 用 于 蒙特 卡 罗 程 序 的 不 确定 度量 化 。 


核 素 及 反应 
道 
Nuclides and 


recations 


表 3 基于 多 群 协 方差 和 连续 能 量 协 方差 的 相对 不 确定 度 计算 结果 对 比 
Table3 Comparison of relative uncertainty results based on multi-group covariances and continuous- 
energy covariances 


Godiva 基准 题 U-235 
Godiva Benchmark U-235 
连续 能 量 协 
多 群 协 方差 ”方差 计算 结 
计算 结果 果 
multi-group continus 
covariance energy 
calculation covariance 
results calculation 
results 
2.900E-03 2.920E-03 
3.990E-03 3.995E-03 
2.750E-03 2.790E-03 
7.880E-03 7.880E-03 


本 文 针对 评价 核 数据 库 


介 协 方差 数据 的 处 理 ， 


核 数据 处 理 程序 NECP-Atlas 中 开发 了 
自主 化 ， 实 现 了 平均 裂变 ! 


给 出 的 各 类 型 的 协 方差 数据 ， 分 别 而 


子 数 、 共 振 参数 、 中 子 截面 、 


产生 对 应 的 多 条 


库 的 计算 


结 


应 根据 ， 


子 学 输 运 计算 程序 
应 选择 与 多 群 数据 库 版 本 一 致 、 能 条 


相对 偏差 /% 


realative 


bias/96 


0.69 
0.01 
1.45 
0.00 


Jezebel 基准 题 Pu-239 
Jezebel Benchmark Pu-239 


群 协 方差 数据 库 进 行 比较 以 及 基准 题 的 检验 , 验证 了 NECP-Atlas fi 


果 对 比 , 证 明了 连续 能 量 协 方差 数据 
的 特点 ， 制 作 并 使 


库 的 正确 性 。 后续 
相应 的 协 方差 数据 库 。 基 于 确 


连续 能 量 协 
多 群 协 方差 ” 方差 计算 结 
核 素 及 反应 计算 结果 ES iu 
i . j 相对 偏差 /% 
道 multi-group continus 
realative 
Nuclides and covariance energy i 
bias/% 
recations calculation covariance 
results calculation 
results 
O, 1.550E-02 1.544E-02 -0.39 
O` hel 8.700E-03 8.670E-03 -0.34 
O; 9.050E-03 9.060E-03 0.11 
V 3.180E-03 3.170E-03 -0.31 
究 了 相应 的 应 用 协 方差 数据 计算 方法 ， 在 
协 方差 计算 模块 covar_calc， 完 成 了 协 方差 数据 计算 模块 的 国产 化 、 
次 级 粒子 角度 分 布 、 裂 变 谱 、 中 子 活化 截面 等 评 
f 协 方差 数据 库 或 连续 能 量 协 方差 数据 库 。 通 过 与 NJOY21 制作 的 多 
制作 的 多 群 协 方差 数据 库 的 正确 性 , 计算 


精度 与 NJOY21 相当 ; 通过 将 NECP-Atlas 制作 的 连续 能 量 协 方差 数据 库 计 算 的 数 


=] 


BARSE EVI E 
究 人 员 开展 不 确定 度量 化 相关 工作 时 ， 


定论 程序 的 不 确 


定 度量 化 ， 


结构 一 致 的 多 群 协 方差 数据 库 


可 直接 采 


作者 贡献 声明 OUR: 程序 开发 、 验 证 计算 、 结 果 分 析 、 起 草 文章 ; 祖 铁军 : 程序 开发 、 理 
及 文章 整体 设计 及 修改 ， 曹 良 志 : 到 


版 本 一 致 的 连续 能 量 协 方差 数据 库 。 


neS 


RRF: 


ER 
Js 


里 论 指导 。 
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定 度量 化 ， 


H 


论 指导 ， 以 
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